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SUMMARY 

An  exact  solution  is  found  for  the  temperature  distribution  in  a 
rectangular  plate  of  uniform  thickness  under  varying  ambient  temperature. 
For  the  special  case  of  temperature  variation  in  the  form  of  a  sinusoidal 
pulse,  a  numerical  comparison  is  made  of  the  exact  solution  and  an  ap¬ 
proximate  solution  obtained  by  a  modification  of  the  solution  for  con¬ 
stant  ambient  temperature.  The  range  of  validity  of  this  approximate 
solution  is  considered  in  terms  of  a  "rise  time." 
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PHASE  REPORT 

HEAT  PLOW  IN  RECTANGULAR  PLATE  OF  UNIFORM 
THICKNESS  UNDER  VARYING  AMBIENT  TEMPERATURE 

WEPTASK  NO.  RAViilJ002/202l/R008-03-001 
Problem  No.  3  (AV-liij.Ol6) 


Ref:  (a)  Adams,  Edwin  P.  and  Hippisley,  R.  L.  -  Smithsonian  Mathemati¬ 
cal  Formulae  and  Tables  of  Elliptic  Functions,  Washington  D.C., 
1939 

(b)  Report  No.  NADC-EL-6l[i3,  "Phase  Report,  Development  of  Analyt¬ 
ical  Design  of  Heat  Exchangers  for  Airborne  Electronic  Equip¬ 
ment" 

INTRODUCTION 


The  problem  of  calculating  the  cooling  effect  of  fins  under  steady 
state  conditions  is  relatively  elementary,  in  that  it  involves  the  solu¬ 
tion  of  an  ordinary  differential  equation.  However,  under  non-steady 
state  conditions,  as  when  the  ambient  temperature  is  varying,  a  partial 
differential  equation  is  involved.  The  solution  to  this  equation  must 
satisfy  a  specified  initial  temperature  distribution  as  well  as  the 
boundary  conditions. 

The  rate  of  heat  flow  into  the  base  of  the  fin  in  the  steady  state 
may  be  written  as 

Q  =  (T^  -  T3^)F(P,G) 

where 

=  temperature  at  base  of  fin 

T^  =  ambient  temperature 

F(F,G)  =  a  function  of  the  physical  and  geometrical  parameters  of 
the  fin. 


If  now  the  ambient  temperature  is  a  variable,  one  is  tempted  to  avoid 
the  difficulties  associated  with  the  solution  of  a  partial  differential 
equation  and  to  use  the  above  equation,  derived  under  the  assumption  of 
steady  state  conditions,  by  substituting  the  variable  ambient  temperature 
Tg^(t)  for  the  constant  ambient  temperature  T^,  i.e.,  by  writing 


Q(t) 


Tjt) 


F(P,G) 


and 


Q(t)  =  F(P,G)/  |t^  -  T^(t) 


dT 


This  approximation  will  be  referred  to  as  "quasi-static." 
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There  is  no  doubt  that  thi '  procedure  furnishes  a  good  approximation 
if  Ta(t)  varies  slowly  enough;  but  to  justify  its  use  the  expression 
"slowly  enough"  must  be  made  more  quantitative. 

In  this  report  the  general  solution  of  the  proper  partial  differen¬ 
tial  equation  for  a  fin  with  fixed  base  temperature  and  variable  ambient 
temperature  is  derived.  This  is  then  specialized  to  the  case  where  the 
ambient  temperature  varies  as  a  step  function  or  as  a  sinusoidal  pulse. 
For  this  latter  case  the  exact  solution  is  compared  with  the  solution 
obtained  by  the  approximate  method  discussed  above,  and  the  difference 
is  discussed  in  terms  of  a  "rise  time"  obtained  from  the  response  to  a 
step  function. 

For  simplicity  the  fin  was  assumed  to  have  a  uniform  thickness 
(smallest  dimension),  but  the  general  conclusions  are  applicable  to  the 
case  of  a  tapered  thickness. 

THEORY 


Differential  Equation 


/) 


/ 


■1 


//  unit  width 
/  / 


XX 


d 

A - y 

T 


X  =  0  dx  X  ^  -  L 

Assume  heat  flov/  in  x- direction  only.  Then, 


net  rate  of  heat  input  to  element 


rate  of  heat  transfer  tc  ambient  it 


medium  =  2h|^T  -  T^(t) 


dx 


ST 

rate  of  heat  accumulation  in  element  =  cp6  dx 


^^2  0 

k  Bt 


Let  a 


3  _  2hL^ 


6  k 


3  =  =  T„  +  f(t);  then 


2 
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(1) 

(2) 

(3) 


Equation  (3)  expresses  the  assumption  that  the  free  edge  of  the  fin 
is  insulated,  so  that  no  heat  flows  through  it. 

When  t  =  0  the  temperature  distribution  assumed  is  to  be 


T(0,x)  =  Ta  +  (Tt,  -  Ta) 


cosh  a(l  -  x/L) 
cosh  a 


This  is  the  steady  state  temperature  distribution  with  ambient  tempera¬ 
ture  T^.  (This  temperature  di striluti cn  is  easily  derived  from  the 
steady  state  heat  flow  equation.) 

General  Solution  by  Means  of  Laplace  Transform 

Let  6 (s,x)  =  oZC  [ T(t,x)] 

F(s)  =/-{f(t)} 

where 

X [  }  =  Laplace  transform  of  {  ] . 

Then 

X  =  se(s,x)  -  T(0,x) 


and  with  the  use  of  (it)  the  transform  of  equation  (l)  becomes 

(a  +0s)9(s,x)  =  -  —(a  +  g s )  -  F( s )  -  g ( Tb  -  Ta )-^  —  ~ 

d(x/L)  s  cosh  a 


The  general  solution  of  this  ordinary  differential  equation  is 
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6(s,x)  =  A  cosh^a^  +  3s  ^  +  B  sinh^,  +  3s  ^ 

^  a^F(s)  ^  £a  ^  _  ^^^cosh  a(l  -  s/L) 

+  3s  s  s  cosh  a 

The  transformed  boundary  conditions  become 


9(s,0)  =  Tb/s 

(2) 

o 

II 

W 

CD 

(3) 

Substitution  of  (5)  in  (2)'  gives 

_  a"F(s) 

“  +  3s 

Substitution  of  (5)  in  (3)'  gives 

B  =  -  A  tanh^OL^  ■'■3  s 

Thus 


0(s,x) 


a^F(s) 
+  3s 


>  Q  ^  Q  ^  O  JV 

cosh^'a  +  PS  -  -  tanh^^cT  +  3s  sinh.ya  +  3s  - 
1;  L 


.  F(s)  ^ 
+  3s 


la  1 

~  ^  -(Tb  -  Ta) 

s  s 


cosh  a(l  -  x/L) 
cosh  a 


or 


0(s,x) 


'F(s) 


+  3s 


cosh^.  (£"  +  3s  (l  -  x/L) 
cosh  +  3s 


"  -  Ta  MTb  -  Ta 


^cosh  a(l  -  x/L) 
cosh  a 


Now 

T(t,x)  =y;“'[e(s,x)] 
where 

]  =  inverse  Laplace  transform  of  [  }  . 
The  first  term  may  be  written  as  F(s)G(s);  and  if 


-  h  - 
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T(t,x)  =  -  sin  m“  ~(l  -  e 

n  +  a"  )  2  L 


2TT^  ,  2 


n  X,  -  (nf"  +  a^)t/0, 
I - (1  -  e  4 


"  -  4) 


cosh  a(l  -  x/l) 
cosh  a 


The  rate  at  which  heat  enters  the  base  of  the  fin  (per  unit  width) 
from  the  body  being  cooled  is 


,  ST 
Q(t)  =  -6k — 


Sxlx  =  0 


2a  1  . 

''  1,3,-  ”  X  “ 


)t/3 


+  k6^(T^  -  T^)  tanh  a 


This  can  be  written  in  a  normalized  form  as 

A  ^  Q(t)  -  Q(0)  ^  2a  1  -  ^  )Vp 

4  (t)  =  T-— - r—  =  1  -  -  /  S - re 

Q(“)  -  Q(0)  tanh  a  m  ^  (p) 

where  use  has  been  .'de  of  the  relation 


UJ 

1,3,- "T 


1  tanh  a 


2  j.  2  /-) 

m  -I-  +  a  2  a 


(This  equation  follows  from  equation  6. [[95-1  in  reference  (a);  but  it  can 
also  be  obtained  from  the  condition  that  if  AT  =  4  “  4’  Q(“)  =  0.) 

2.  (t)  is  zero  for  t  =  0  and  approaches  unity  asymptotically  as  t  “j 
also,  except  for  values  of  t/3  very  close  to  zero,  it  can  be  expressed 
as  a  simple  exponential  rise 


2(t)  1  -  - 


-  )t/3 


4-  U  -L  2 

tanh  a  —  +  a 


This  equation  is  plotted  in  figure  1  for  several  values  of  a. 
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Thus  the  response  of  the  cooling  fin  to  a  step  change  in  ambient 
temperature  can  be  simply  characterized  by  a  rise  time,  i.e.,  the  time 
required  for^(t)  to  reach  some  fraction,  f,  of  its  final  value,  unity. 
The  rise  time  as  calculated  from  the  above  equation  is  given  in  figure  2 
for  f  =  80  and  90  percent. 

Response  to  Sinusoidal  Pulse 


"Exact"  Solution 


Let 


f(t)  =  { 


AT  sin  -jr 
0  ^ 


0  <  t  <  D 
t  >  D 
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T(t,x)  =  ^ATa^r 


CO 

E 

1,3,- 


1  /  2  TT^  ,  2  \2  .  2 

.m|(in  +  a  j  +  r 


n  X 
sin  ,1712  L 


X 


2^  2 
m  +  a 


TTt 


sin  -  cos  3"  e 


2  n^  ^  2  \Trt 
TTt  T  ^  D 


+  T  +  (T  -  T  )  a(l  -  x/L) 

a  ^  b  a'^  cosh  a 


t  <  D 


jjATa  r 


.  °°  r 

L  ^ 

1,3,- 


+  o'  )' 


,2  ^2  ^  J3 


n  X 
sin  m2  l 


X 


1/  2n^  ,  2  N^(t  -  D)  1,  _2  .nt 

-(m  -j-  +  a  ) - p -  ,  -  T 


cosh  a(l  -  x/L) 


cosh  a 


t  >  D 


where 


r' 


np 

D 


From  Q(t) 


0 


-  9  - 


NADC-EL-6199 


k6Q(t)  - 


-  2ATa‘'r 


00 

E 


x:x 


1  f  3  TT^  ,  2  \  TT  t* 

nt  ^  T  “  V 

COS  -^  +  e  ^ 


+  a(T^  -  Tg^)  tanh 


t  <  D 


,3  y 


1,3,- 


/  2  "  ,  2  \2  , 
(m  -]—  +  a  }  + 


+  a(T|^  -  Tg^)  tanh  a 


t  >  D 


Quasi-Static  Approximation 

From  the  temperature  distribution  for  constant  ambient  temperature, 
as  given  by  equation  (h),  the  rate  of  heat  flow  into  base  of  the  fin 
under  this  condition  is 

^  ■  -  “six  -  0  '  “E<''b  -  ‘a>  “ 


In  accordance  with  the  discussion  in  the  Introduction  an  approximate 
solution  for  variable  ambient  temperature  can  be  written  by  replacing 


by  T^(t).  In  the  case  where 


T,(t) 


m  TTt 

Ta  .  AT  sin  ^ 


for  t  S  D 
for  t  >  D 


the  approximate  solution  becomes 
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i^Q(t)  =  -  atanha  AT  sin  ^  +  (T^  -  a  tanh  a 


k6 


(Tb  -  Tg)  a  tanh  a 


t  <  'D 
t  >  D 


(13) 


DISCUSSION 

The  change  in  cooling  rate  produced  by  a  change  in  ambient  tempera¬ 
ture  in  the  form  of  a  sinusoidal  pulse  was  calculated  from  equations  (12) 
and  (13)  for  three  fins  of  the  dimensions  given  below. 

Fin  1  2 

L  (ft)  0.0625  0.110 

6  (ft)  0.0015  0.0021 

The  assumed  physical  constants  are: 
h  =  1  btu/(hr  ft^  "F) 
k  =  100  btu/(hr  ft  "F) 
p  =  168  Ib/ft^ 

CV  =0.23  btu/lb  °F 

and  the  temperatures  involved  are 
Tb  =  150^  F 
=  95'  F 
AT  =  50"  F 


3 

0.19ii 

0.00327 


The  same  three  values  of  pulse  duration  vrere  chosen  for  each  fin. 

D  (hr)  =  0.22,  0.033&,  and  0.00375. 

The  derived  constants  are : 


Fin  1 


2  3 


a  0.228  0.3ii0 

p  (hr)  0.00151  0.00i|68 

(90,^)  t^  (hr)  0.00125  0.00378 


0.i;80 

O.Oli;58 

0.0112 


The  results  of  these  calculations  are  shown  in  figures  3)  and  5. 
The  ratio  of  the  instantaneous  rate  of  heat  flow  into  the  base  of  the 
fin  to  the  steady  state  rate  is  plotted  versus  the  ratio  of  time  to  the 
endurance  of  the  pulse. 
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The  ratio  of  the  pulse  duration  to  the  rise  time  of  the  fin  (D/tj.) 
is  indicated  on  each  curve.  This  ratio  is  a  measure  of  the  "slowness" 
of  the  ambient  temperature  change,  and,  as  seen  from  the  figures  and  as 
is  intuitively  obvious,  the  exact  solution  approaches  the  quasi-static 
approximation  for  large  values  of  D/t^,. 

The  ambient  temperature  time-history  given  in  figure  6  of  reference 
(b)  may  be  approximated  by  a  sinusoidal  pulse  of  duration  D  =  0.22  hour. 
Thus  curve  A  of  figure  3  gives  the  exact  solution  for  a  problem  analogous 
to  the  one  considered  in  reference  (b).  In  thj-s  case  use  of  the  approx¬ 
imate  solution  is  very  satisfactory,  and  the  small  errors  it  does  produce 
are  nearly  all  on  the  conservative  side. 
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FIGURE  U  -  Fin  2  Heat  Flow  Rate 
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FIGURE  5  -  Fin  3  Heat  Flew  Rat 
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AFPEN'DIX  A 
LIST  OF  SYMBOLS 


L 

6 

P 

Cp 

k 

h 

t 

T 

AT 

Tb 

T 

a 

T^(t)  =  T^  +  f(t) 
f(t) 

D 

a  s  v  2hL^/6k 
P  =  CgpLVk 
r  =  rrS/D 


length  of  fin 

thickness  (smallest  dimension,  sometimes  called 
"width"  in  heat  transfer  literature) 

density- 

specific  heat 

themal  conductivity 

heat  transfer  coefficient 

time 

rise  time 

temperature 

temperature  increment 

base  temperature 

constant  ambient  temperature 

varying  ambient  temperature 

variable  component  of  ambient  temperature 

duration  of  sinusoidal  pulse 

dimensionless  constant 

constant  (dimension  of  time) 

dimensionless  constant 
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APPENDIX  B 
SUMMATION  OF  SERIES 


to 


T7 1 

In  equations  (ll)  and  (12)  the  coefficient  of  sin  is  proporticna 


f  L 


2^7“  2 


.a  72  ^  ^2 


For  large  values  of  r  the  calculation  of  many  terms  in  this  series 
is  required  for  accuracy:  however,  this  series  can  be  summed  exactly. 
It  is  convenient  to  evaluate  simultaneously  the  series 


appearing  in  the  coefficient  of  cos  ^ 


Let 


A  = 


CC 


2'^  2 

T  ^  “ 


^2 

/  2  '  ,  2  \2  ,  .2 
_ (m  —  +  aj  +r 


B  = 


y 


A  +  iB  -  “ 


m 


T 


+  cw=*  +  ir 


(m^^  +  )  +  r^ 


i  r 


m  +  a 


ir 
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X 


y  =  \/ 


r 

/  „ 


a 


2 


yd  ■*■ 


1  u  sinh  2u  +  V  sin  2v' 
2v\/ <x  id  cosh  2u  +  cos  2v 


CD 

1  ■'  V  sinh  2u  -  u  sin  2v 

!■  I  M  c  ■  I  •  ■■  ■  I  ■  I  II  ■  I  . . .  nm  ■■  III, HI  ■■  lll» 

2rv^^  +  cosh  2u  +  cos  2v 

1,3,--  ^ 

where 


u  = 
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